*********************************	
/* Table 2 - Base Models       */
*********************************

* Note it is necessary to install the firthlogit add-on
* findit firthlogit

eststo clear

set more off

/* Model 1 - no selection all cases */
logit absout i.ccp##c.newpow bilat ww2 contig totalmids, cluster(ndyadid) nolog

	estimates store model2

/* Model 2 - no selection drop ww2 */
logit absout i.ccp##c.newpow bilat contig totalmids if ww2 == 0, cluster(ndyadid) nolog

	estimates store model2

/* Model 3 - firth logit all cases */
firthlogit absout i.ccp##c.newpow bilat ww2 contig totalmids

	estimates store model3

/* Model 4 - firth logit drop ww2 */
firthlogit absout i.ccp##c.newpow bilat contig totalmids if ww2 == 0

	estimates store model4

* Table 2 Output
estout model1 model2 model3 model4, ///
	cells(b(star fmt(%9.3f)) se(par fmt(%9.3f))) style(tab) ///
	stats(N, fmt(a1 %9.3f %9.3f)), using "absout_base.txt"







	
